Budding and Domain Shape Transformations in 
Mixed Lipid Films and Bilayer Membranes 
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We study the stability and shapes of domains with spontaneous curvature in fluid films and mem- 
branes, embedded in a surrounding membrane with zero spontaneous curvature. These domains 
can result from the inclusion of an impurity in a fluid membrane, or from phase separation within 
the membrane. We show that for small but finite line and surface tensions and for finite sponta- 
neous curvatures, an equilibrium phase of protruding circular domains is obtained at low impurity 
concentrations. At higher concentrations, we predict a transition from circular domains, or caplets, 
to stripes. In both cases, we calculate the shapes of these domains within the Monge representa- 
tion for the membrane shape. With increasing line tension, we show numerically that there is a 
budding transformation from stable protruding circular domains to spherical buds. We calculate 
the full phase diagram, and demonstrate a two triple points, of respectively bud-flat-caplet and 
flat-stripe-caplet coexistence. 

PACS numbers: 68.18.Jk Phase transitions 



I. INTRODUCTION 



Fluid membranes occur in a wide variety of physi- 
cal, chemical, and biological systems Q, H 13, 3 Ex- 
amples include surfactant films, uni-lamellar and multi- 
lamellar vesicles, and lipid bilayer membranes, such as 
those in biological cells. Increased attention has been 
paid to the properties of multi-component bilayer mem- 
branes. There are at least two important reasons for 
this. On the one hand, biological cell membranes natu- 
rally involve mixtures of several different lipid and pro- 
tein components [1 • Processes such as budding, shape 
and textural transformations, and raft formation are be- 
lieved to involve local inhomogeneities in these mixtures 
H S 0, S ©, G3 [H El El M El El E3- On the other 
hand, lipid mixtures have been shown to play an impor- 
tant role in the formation of stable vesicles for a variety of 
potential biomedical applications such as controlled gene 
and drug delivery |l8l [l9j • 

The additional internal degrees of freedom that accom- 
pany the presence of two or more components within a 
fluid membrane can lead to a rich set of different mem- 
brane properties. The understanding of the effects of 
composition variations on membrane properties is im- 
portant, as it may shed light on the behavior of cell 
membranes, and also enable rational control of synthetic 
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membrane structure and function. The compositional 
degrees of freedom in a multi-component membrane can 
dramatically influence both its morphology and phase be- 
havior. For instance, phase separation within the mem- 
brane can occur, analogous to phase separation of ordi- 
nary fluids. However, such two-dimensional phase sep- 
aration within a membrane is expected to be closely 
linked with the three-dimensional shape of the mem- 
brane. This is because, at a microscopic level, the cur- 
vature properties, and especially the spontaneous curva- 
ture, are largely dictated by the structure of the con- 
stituent molecules 20] . Thus, regions of different com- 
position will often have different curvature properties. In 
such cases, the formation of a domain of one phase within 
a matrix of another phase can lead to a localized deforma- 
tion of the membrane. This deformation can be enhanced 
by finite line tension between domains, which when suffi- 
ciently large can drive budding (the formation and subse- 
quent separation of a small vesicle). Experimental stud- 
ies of model two and three component systems have ex- 
amined domain formation HHHHE3E1E1E1E1 
and its role in the process of membrane shape transfor- 
mations 0, 0| and budding U Q • 

Many theoretical studies have examined the interplay 
between internal degrees of freedom and membrane shape 
in two-com pon ent fluid membranes l2ll I22I I23I Izil EM 
ElllllllllllMIMIMIlllMlIl. In theoretical 
studies of phase-separated bilayer vesicles, it was found 
that budding can occur in the limit of large line ten- 
sion between the two phases |28|, |2!j, |33j. However, 
when surface tension is relevant, phase separation can 
lead to stable, modulated phases of flat films and vesicles 
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mmiHEllillUHIil. Indeed ' equilibrium mod- 
ulated morphologies have been observed in several re- 
cent studies of phase-separated multi-component bilayer 
vesicles j^, ■ O n the other hand, the deformation of 
homogeneously mixed two-component fluid membranes 
may induce in-plane phase separation 30J . 

For a mixed monolayer, the deformation of the film 
due to phase separation of two components with differ- 
ing spontaneous curvature is an immediate consequence 
of the different molecular architectures of the two com- 
ponents. For a bilayer membrane, however, the coupling 
of composition to curvature is somewhat more subtle be- 
cause of symmetry considerations. Phase separation can 
take two different forms within a bilayer membrane. As 
shown in Fig. ^a), a symmetric domain of one com- 
ponent can develop within a matrix of another phase. 
In Fig. Q[b) , an asymmetric domain is shown, in which 
the two halves of the bilayer have different compositions. 
As illustrated in the figure, these effects can influence 
the shapes and phase behavior of membranes in differ- 
ent ways. In the symmetric case [Fig. [Tfa)], the domain 
can remain flat |28| although the domain can be more 
or less rigid than the surrounding membrane, a situation 
reminiscent of rafts in biological membranes. However, 
composition inhomogeneities can give rise to deforma- 
tion of a bilayer membrane if there is an asymmetric dis- 
tribution of the lipid constituents across the bilayer, as 
shown in Fig. ^b) |36| . Such lipid asymmetry has been 
clearly demonstrated experimentally in biological mem- 
branes [37L l38| . This broken symmetry of the bilayer can 
either arise spontaneously [31| or as a result of different 
environments on the inside and outside of a vesicle. In 
the latter case, we must also consider the possibility of 
a non-zero surface tension due to an osmotic pressure 
difference between the inside and outside. 



A B 




FIG. 1: (a) Flat symmetric phase-separated membrane, and 
(b) curved, asymmetric phase-separated membrane. 

In this paper, we examine the stability and equilibrium 




FIG. 2: Schematic representation of phase-separated circular 
(top) and stripe (middle) domains and a bud (bottom) of the 
B component, with induced deformations. 



shapes of domains in asymptotically flat fluid films (or gi- 
ant vesicles with dimensions much larger than the domain 
sizes). We focus on domains of constituents with finite 
effective spontaneous curvature, embedded in a matrix 
of membrane material with zero spontaneous curvature. 
These domains can either be due to phase separation 
of two or more components in a mixed membrane, or 
to the inclusion of an "impurity" in the film (such as a 
membrane protein or a surface adsorbed macromolccule) 
[3^ |. We show that for membranes under tension, sta- 
ble protruding circular domains (we call these caplets, as 
in Refs. |32j,|35|) similar to those observed for mixtures 
of lecithin and phosphatitic acid [3^| can occur at low 
concentrations of the minority component or impurity. 
We illustrate the structure of such circular caplet do- 
mains in Fig. |3 We also demonstrate the possibility of a 
phase transition from a mesophase of circular caplets to 
a mesophase of protruding stripe domains either with in- 
creasing area fraction of the minority component or with 
increasing tension. We illustrate these stripe domains in 
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Fig.H 

In both cases, the domains are shown to be stable with 
respect to complete phase separation in a flat membrane. 
In Section^we present our model and briefly discuss our 
analysis in the Monge gauge treatment. In Section ITTT1 
we provide the resulting Euler-Lagrange equations and 
boundary conditions used in our Monge gauge calcula- 
tions. In Section IIVI we present Monge limit results in 
the dilute and concentrated regimes, and discuss the ex- 
tent to which the limiting approximations of Ref. |32l | 
apply. In Section^we present the results of a numerical 
calculation of the budding transition beyond the Monge 
approximation, and we conclude in Section IVII with a 
summary and discussion. 

II. MODEL (MONGE GAUGE) 

We consider a single idealized two-component mem- 
brane consisting of incompatible amphiphilcs A and B. 
Furthermore, we consider asymmetric bilayers, in which 
phase separation occurs in only one leaflet of the bilayer, 
as shown in Fig. ^b). This applies, for instance, to the 
case of phase separation in an asymmetric bilayer vesicle 
or biomembrane. Such domains can also be induced by 
the adsorption of macroions on the inside or outside of 
a bilayer vesicle, as in the experiments of Refs. fft liof. 
We assume that the binary liquid mixture is far from its 
critical point, implying that there are sharp interfaces be- 
tween the domains of A and B, and that the membrane 
consists of a fluid minority phase B surrounded by an 
asymptotically flat fluid majority phase A. 

We investigate the properties of two prototypical do- 
main morphologies (see Fig- EJ) , which are idealizations 
of actual structures that may exhibit disordered textures 
of polydisperse, irregular domains: 

i) Quasi-one-dimensional domains of B (stripes), with 
an impurity phase of width 2xq embedded inside a 
majority stripe with total with 2Xq. The stripe 
length (system dimension) is L, and the projected 
area of an element (minority plus majority repeat 
unit) is Sa.± + Sb,±_ = 2XqL. 

ii) Monodisperse circular domains of B (caplets). For 
caplets we take a circular Wigner-Seitz cells with 
impurity radius tq and total radius Rq. The circu- 
lar Wigner-Seitz cell approximation leads to a lower 
bound for the free energy of a less symmetric struc- 



ture, such as a hexagonal array. The projected area 
of an element is Sa.± + Sb _l = vri?Q. 

Below we will typically scale these dimensions by the 
spontaneous curvature Co, thus defining po = tqCq, Rq = 
Rqcq, £o = £oCch So = XqCq. 

For given material parameters, we obtain minimum 
energy height profiles h(r) (within a linearized Monge 
gauge description), where r is the position in the plane, 
and minimize the resulting free energy over the domain 



(a) (b) 




FIG. 3: A sketch of the stripe and cap let morphologies. 
Fig. (a) shows a linear array of B-phase stripes of width 2xo 
(black) in a matrix of phase A (grey). Fig. (b) shows a hexag- 
onal array of B-phase circular caplets of radius vq (black) in 
a matrix of phase A (grey). The dashed lines in each fig- 
ure show the appropriate Wigner-Seitz cell boundaries, while 
the membrane displacement profile of an appropriate unit cell 
cross section is shown above each figure. 

dimensions. We consider the background A phase to be 
flat and restrict our attention to the simplified case of 
equal bending moduli ka = = k. 

For phase separated membranes under tension, the free 
energy of a single domain (stripe or caplet) has the form 

G =F l +F s +F b (2.1) 

where Fi, F s , and Ff, are the energy contributions from 
the interfacial (line) tension between the two phases, the 
membrane frame (surface) tension, and the bending elas- 
ticity. To lowest order in the Monge representation (see 
Appendix^), valid for nearly flat membranes, these con- 
tributions relative to a flat reference state are given as 
follows. 



The bending energy is 



Sb 



(V 2 h) 2 -2c V 2 h + c 2 (Vh)' 



dS A 



(2.2) 



where dS± is the (flat) area measure in the plane, Si denotes the projected area occupied by phase i, h(r) is the 
height profile above some flat reference plane, V is a two-dimensional gradient operator, k is the bending modulus, 
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and Co is the spontaneous curvature of the B phase. The line tension energy is 

r ( X2nr N (caplets) 

Fi = Xd> dl=l (2.3) 
Jas B { X2LN (stripes) 

where N is the number of domains, 8Sb denotes the boundary of the impurity region, and A is the line tension. The 
work done by the membrane in deforming relative to the flat state is 



(2.4) 



where a is the frame tension, and the Lagrange multiplier 
(or exchange chemical potential per unit area) (i controls 
the amount of impurity phase. We assume a fixed area 
per molecule in the membrane, and thus ignore stretching 
energies. 

It will be convenient to define a renormalized chemical 
potential p,, by 



(2.5) 



where kcq arose in Eq. I|2.5|l from keeping all terms to 
second order in the height field |35|. This renormaliza- 
tion, neglected in |32| . is typically very small: kcq is of 
order fc^T divided by a curvature radius squared of order 
(lOOnm) 2 , while a ~ fc^T divided by a length squared 
of order nm 2 . In fact, we will find that this renormal- 
ization does not change the phase boundaries if o nly one 
phase has a non-zero spontaneous curvature; see [35j for 
a different example. 

Stationarity of Eq. (|2.1|) with respect to height vari- 
ations leads to Euler-Lagrange (E-L) equations for the 
height profiles in regions A and B. The detailed deriva- 
tion of the these equations are presented in Appendix lAl 
To calculate the phase behavior we must minimize the 
free energy G of an array of domains, 



where 



G(ji ) <T,\)=NG (p ) <T,\), 

Sf r 



N 



Sa,± + Sb . 



(2.6) 



(2.7) 



is the number of domains and Sf r is the frame area, which 
without loss of generality we consider to be fixed. The 
free energy must also be minimized over the dimensions 
of the A and B domains, Xq and Xq for the stripes and 
ro and Rq for the caplets: 



dG 
dx 
dG 
dr 



0, 



0, 



dG 
dXo~ 
dG 



= 0, (stripes) (2.8a) 
= 0, (caplets). (2.8b) 



Physically, minimization over the inner length (xo or ro) 
is equivalent to allowing exchange of B for A species, 
while minimization over the outer length (Xq or Rq) is 



equivalent to allowing more or less total area into the 
system, doing work against the frame tension. 

Phase coexistence is found by equating chemical po- 
tentials, 



/^stripe 



/^caplets- 



(2.9) 



In the dilute limit, where all the impurity is in a single 
domain, we must minimize the free energy per impurity, 



A) 



Gg(fi,a, A) 



(2.10) 



over the impurity domain size. In using the projected 
area of the B phase we have kept terms consistent with 
the Monge gauge calculation (and hence ignored a term 
of order h 2 in the denominator). 

Following 1^3, we will use hatted variables for di- 
mensionless quantities A, a, ji (introduced below and Ap- 
pendix^) scaled by n and appropriate powers of cq, and 
unadorned variables for physical quantities. This differs 
from the more recent work of [35| in which dimensional 
quantities were hatted and dimensionless quantities were 
not hatted. 

For simplicity, we have also assumed that all bending 
moduli are the same in the different phases. Different 
mean curvature moduli k will obviously shift the phase 
boundary in favor of the phase with a lower curvature 
modulus, while a difference in Gaussian curvature moduli 
would lead to a shape-dependent line tension that will 
shift the phase boundaries [33|. 



III. EULER-LAGRANGE EQUATIONS AND 
BOUNDARY CONDITIONS 



Variation of Gq([i, a, A) with respect to h(r) gives the 
Euler-Lagrange equations for the height profile. To cal- 
culate the free energy we only need the slope r/ — dh/dx 
(stripe) or r\ — dh/dr (caplet). As shown in Appendix lAl 
the slopes in regions A and B (t = A, B) satisfy 



de'^ ri = 



d 2 rj 
d^ 2 



1 drj 
pdp 



(stripes) (3.1a) 
(caplets) (3.1b) 
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where we have introduced dimensionless variables, 



A. Dilute Limit (a » fj, ~ 0) 





= c x 


(3.2a) 




= c r 


(3.2b) 




= u + (i 


(3.2c) 


ci 


= a 


(3. 2d) 


a 


a 

2 

KC Q 


(3.2e) 


p, 


= A + l 


(3.2f) 


A 


A 

KC ' 


(3-2g) 



Note that the renormalization of p due to kcq is identical 
in both stripe and caplet phases, and will thus play no 
role in determining coexistence, i.e., in Eq. I|2.9f) . 
Symmetry of r/ at the domain centers implies 



77(0) = r/'(0) = 0. 



(3.3) 



where primes denote derivatives with respect to either p 
or £. Another condition may be supplied by requiring 
a well-defined energy. A discontinuity in rj would imply 
singular curvature at the boundary, leading to a non- 
physical curvature energy in the region of the interface 
ds B between A and B phases. Hence we must have 



0. 



(3.4) 



The final boundary conditions follow from the variation 
of the height profile, and are derived in Appendix ^] In 
the dilute limit, 



[cb - c a \ 9sb = c 
(T77a(oo) — kc a (co) = 
ca(oo) = 0, 

where c(r) is the curvature: 



1 d 



(pv) 



(torque) 

(normal force) 
(torque) 

(stripes) 
(caplets). 



(3.5) 

(3.6) 
(3.7) 



(3.8) 



, pdp 

Eqs. 13.5|l and l|3.7|l are torque balances at the A— B inter- 
face |4l[ and at the system boundary, and Eqs. (|3.6[) are 
the consequence of zero vertical force at the A — B inter- 
face. For the concentrated limit, the last two equations 
should be replaced by the equivalent result for periodic 
boundary conditions: 



Was, 



= 



(concentrated regime). (3.9) 



IV. MONGE LIMIT: RESULTS 



The general solutions to the Euler-Lagrange equations 
are given in Appendix [Bj Before discussing the phase 
behavior in the general (concentrated) case, we discuss 
the dilute limit for large frame tensions, where several 
approximations simplify the results. 



Since the height profiles decay in the A phase with a 
characteristic length ^/k/ct, the dilute limit is valid pro- 
vided that the domains are separated by more than this 
length. The general solutions in the Monge gauge, Eqs. 
(|B5llB16|), hav e inverse penetration depths £a = \fb and 
Cb = V" + A- I n the dilute limit and for large frame 
tension we may ignore /j relative to o in the penetration 
depths. Taking the limits Ro —> oo and Xq — > oo, the 
solutions (given in Appendix |B|) reduce to 



fsinh[v / o'C] 
sinhtv^o] 



(£ < Co) 



and 



[exp[V^(£o - 0] (£>&), 



v(p) 

Vo 



where 



Vo = 



= < 



h [Vd-po] 
. Kx [V&po] 



(P < Po) 



(P > Po), 



(stripes) 

(4.1) 

(caplets) (4.2) 




a ( 1 + coth 



(stripes) 



h [Vfrpo] Ki [Vfrpo] (caplets) 



(4.3) 



Substituting these profiles into Eq. (|2.1|) . the free en- 
ergy per impurity of a single domain is (recall that we 
set fi — in the dilute limit) 



Go 











A - 


ho 










2A 


- Vo 


I Po 





(stripes) 
(caplets), 



(4.4) 



where the impurity area Aq is 2Lxo (stripes) or irr, 
(caplets). The free energies can be rewritten as 



Go 



1 + coth v^Cf 



2r - Vfrpoh (Vfrpo) Ki (y&po) 



V^Po 



(stripes) 
(caplets), 

so that the phase diagram in the dilute limit only depends 
on the control parameter t = <7 1 / 2 A |46|. Minimizing 
the free energies over the periods yield the equilibrium 
domain sizes, shown in Fig. 21 The caplet phase is stable 
with respect to stripes for r < 0.2115, while stripes are 
unstable with respect to a flat phase for r > 0.25. 
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FIG. 4: Reduced domain size as a function of r in the di- 
lute limit (</>—» 0,/t —> 0). The solid lines denote the stable 
phase. Caplets are stable for r < 0.2115, and stripes are sta- 
ble for 0.2115 < t < 0.25. For r > 0.25 the stable state is a 
macroscopically phase separated flat phase. 



The stripe domain size diverges as the limit r = 0.25 
is approached, as 



lim fffVa = I ln f^_ 

r^0.25 S 2 V4T-1 



(4.6) 



while the caplet radius grows from zero for small t as 

1/3 



lim p<r 1/2 



T — 









v4r)_ 







(4.7) 



Thus, the domain size diverges as either the line tension 
or the surface tension becomes large. The latter case sim- 
ply corresponds to the fact that for large surface tensions 
there will be complete phase separation in a nearly flat 
two-dimensional system. The former case is somewhat 
more subtle, although, as we shall see, this limit cor- 
responds to budding, since the resulting slope 770 above 
becomes large. 



B. Finite Concentration of the B Phase 

For finite impurity concentration of the two phases, 
the the effective surface tensions of the two phases dif- 
fer by the chemical potential (Lagrange multiplier) fi, 
which cannot strictly speaking be ignored. In this case 
the phase diagram does not depend solely on the dimen- 
sionless quantity r = Act 1 / 2 . This was also noted by 
Gozdz and Gompper, in the case of systems with two 
different spontaneous curvatures (35|. However, we will 
see that the reduction to a phase diagram that depends 
only on r is an excellent approximation in many cases. 
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FIG. 5: Chemical potential as a function of generalized ten- 
sion for different frame tensions a. Note the extrapolation to 
small jX with increasing r. 



For finite area fraction <j> of the B phase, we assume 
a regular array of monodisperse (stable) domains, each 
within cell dimensions determined by the concentration. 
We approximate the hexagonal Wigner-Seitz cell (se e 
Fig. by a set of circular domains of equal area • 
The solutions are given by Eqs. (|B5|) (stripes) and (|B16|) 
(caplets), with corresponding free energies g given by 
Eqs. IjBlOl) and 1B20(I . To calculate the phase diagram 
we minimize the free energies over both domain dimen- 
sions (xo and Xq for stripes, ro and Rq for caplets) for a 
given set of control parameters (line tension A, frame ten- 
sion ct, chemical potential /t), and determine the phase 
boundaries by that chemical potential fi, for which the free 
energies of the caplet and stripe phases are the same. 

Our procedure differs slightly from that of Ref. j^] 
in that, here, the effective overall "tension" is different 
in the two phases. In the impurity phase B the quan- 
tity a + jx acts like a mechanical tension penalizing area 
changes, while (because we have chosen a reference chem- 
ical potential pbA = 0) only ct penalizes area changes in 
the A phase. Hence the inverse penetration lengths (a 
and £b differ. In the strong tension limit ct 3> /t, which 
is exact for vanishing (f> where /t becomes quite small (se e 
Fig. [5J, this approach recovers the results of Ref. |32fl . 
while for smaller ct the height profile of the impurity 
phase contains a contribution from the chemical potential 
as well as the frame tension. In the limit ct 3> £l where 
the penetration depths are the same, the phase diagram 
depends on only the generalized reduced tension r. 

With increasing cf>, the domains grow monotonically 
with ro((f>) ~ 2£o(0) , which corresponds to nearly equal 
mean curvature c ~ Co- Note that a finite preferred 
domain size implies the existence of an equilibrium 
mesophase. With increasing <fi there is a first order tran- 
sition from circular caplets to stripes for small r (Fig. HJJl . 
For small r the transition depends very weakly on r and 
coexistence occurs over a range of concentrations near 
4> = 0.35. This range moves to smaller concentration 
with increasing r, while for 0.2115 < t < 0.25, however, 
the stripe phase becomes stable at all concentrations. 

Figs. El El and [7| show that for large ct the phase dia- 
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: a i/2 

X = AG 



FIG. 6: Phase diagram as a function of B-phase area fraction 
4> and generalized reduced tension r, for different constant 
frame tensions a (A was varied to scan r). The region be- 
tween upper and lower pairs of lines are two-phase regions of 
coexisting stripes and caplets. The inset shows convergence 
for large a near the small r phase boundary. An equilibrium 
flat phase is predicted for r > 0.25. Scans for different a 
converge to the results of Ref. 32] , which corresponds to the 
limit <t — > oo. 



grams converge to a scaling form that depends only on 
r, recovering the phase diagram calculated in Ref. [3^ |. 
For larger a the concentration range of the stable caplet 
phase widens at small r, moving to larger <j>. Recall, 
Eq. l|8.2c[l . that the characteristic inverse decay length 
in the impurity phase is £b = vct + Mj where we find 
jl < along the phase boundary. Hence, for a smaller 
tension a the negative potential decreases the energetic 
cost of buckling against tension, which destabilizes the 
stripe phase. The deviation of the phase boundary near 
t ~ 0.2115 is due to departures from the Monge limit. 

The length scales in the concentrated limit are larger 
than in the dilute limit, for the same r, and converge to 
the dilute limit result when the caplet phase loses stabil- 
ity at t — 0.2115, as expected (Figure|SJ. 



NUMERICAL SOLUTIONS IN THE 
BUDDING REGIME 
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FIG. 8: Scaled domain sizes for caplet and stripe phases along 
the phase coexistence line, for both the concentrated and di- 
lute limit calculations. For each geometry the two dimensions 
(stripe half-period £o and minority domain half width Ho, and 
caplet inner radius po and Wigner-Seitz radius R, all scaled 
by <t) are shown. The variation in domain sizes with tension 
a can be seen to be very small. 



The calculations presented above and in Ref. |32j are 
based on a small-slope approximation, and are only valid 
for small tensions a. In order to both test the validity 
of this approximation as well as to study the possible 
transition from the caplet state above to a budded states, 
we have performed direct numerical minimization of the 
model free energy in Eq. (|2.1|) . beyond the small-slope 
approximation. We present here the equilibrium shapes 



and phase diagram for the stripe phase, the (azimuthally 
symmetric) caplet phase, and for buds. 

We focus on the dilute limit, i.e., for small area fraction 
if). It is expected that phase separation may, under appro- 
priate conditions, lead to the formation of buds (nearly 
spherical domains of one phase connected with a narrow 
neck to another phase that is flat). It has previously 
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been shown theoretically that the line tension between 
two phases alone may result in budding [2^, ■ The 
role of phase separation in budding and fission has also 
been studied in recent experiments |44|- As described 
above, our small-slope approximation used in the previ- 
ous sections is expected to fail for large line tension A. 

In order to address both the validity of the preceding 
analysis, as well as to treat possible budding transitions, 
we describe the membrane shape by a local tangent angle 
9 relative to a horizontal, flat conformation. For a stripe 
domain, the shape can be completely characterized by the 



tangent angle as a function of only one contour length co- 
ordinate s, which is defined along a line perpendicular to 
the stripe. For a caplet or bud domain, a single tangent 
angle as a function of a single coordinate s can also be 
used to completely describe the membrane shape, even 
in the presence of overhangs. This can be done provided 
that the domain is symmetric about an axis perpendicu- 
lar to a horizontal plane. In this case, the contour length 
coordinate is defined along a radial direction. A similar 
approach was used in [3j| . 
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FIG. 9: Data collapse of the free energy per unit area (in units of kcq), domain size (units of Co) and boundary slope as a 
function of r for both stripes and caplets, for various values of the line tension A. 



If the center of the domain is taken as the origin, then 
the membrane shape in three dimensions is given by 

(x(s),y(s),z(s)) = (r(s) cos(</>), r(s) sin(0), z(s)) , (5.1) 

where r(s) = J cos((9) ds is a radial coordinate, z(s) = 
J sin(0) ds, and c6 is the azimuthal angle. In this coordi- 
nate system, the mean curvature is given by 



for caplet and bud domains. For stripes, only the first 
term above is necessary. 

The resulting free energy is evaluated numerically for a 
domain shape defined by a discrete set of 50 to 200 angles 
{0i}i=o,...,AT) which are defined at equally-spaced points 
Si along the contour from the origin to the Wigner-Seitz 
boundary described previously. These values, together 
with an overall scale factor As = Sj+i — Sj determine the 
membrane shape. 
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A discrete approximation to the free energy / per unit 
area of domain is minimized to determine equilibrium 
shape and free energy. Here, we quote results for an area 
fraction of 10%. As shown above within the small-slope 
approximation, the phase boundary between stripes and 
caplcts is very narrow and relatively insensitive to area 
fraction in this range. In Fig. EI we show the normalized 
stripe domain size (£o)i boundary slope (?7o), and free en- 
ergy / vs r for several different combinations of a and A. 
The apparent collapse of the data demonstrates one of 
the conclusions of the previous analysis, namely that the 
results can be well-represented in terms of a single com- 
bination of parameters, r = \^f&, of surface tension and 
line tension. Similar results for caplet domains are shown 
in Fig. |5J Here, too, the results show the dependence on 
the single parameter r. 

These numerical results can also be used to determine 
the phase boundary between the caplet and stripe phases. 
For simplicity, we do this only in the dilute limit, in which 
the phase boundary is rather insensitive to concentration. 
Specifically, we find the crossing of the free energies of the 
caplet and stripe phases (f c and f s ) for 10% area fraction. 
The results are shown in Fig. We note the excellent 
agreement of this phase boundary with that shown in 
Fig. El at least for A < 0.2. For larger A the former 
approximations fail. One consequence of this is that the 
stripe phase vanishes with increasing A: specifically, there 
is a narrow range of values of the line tension for which 
the only stable domains are caplets. This can be seen 
in Fig. 1101 The simple physical reason for the enhanced 
stability of the caplet phase (compared with the small- 
slope results above) in this limit is that the line tension 
cost is lowered in the caplet geometry for large slopes, 
given a fixed domain area: for caplets, the circumference 
to area ratio is reduced with increasing slope. 

Also shown in Fig. ^| is the transition between caplet 
and bud states, as well as a spinodal line beyond which 
no stable caplet structures are possible. No spinodal rep- 
resenting the limit of metastability of buds is found. This 
is an artifact of the use of a bending free energy using 
only leading order terms. Such a spinodal depends on 
the details, and especially the radius a of the neck in 
the budded state. Without further assumptions, this 
length shrinks to of order As. The surface tension re- 
quired to pull-out such a bud is of order the two dimen- 
sional Laplace "pressure" a ~ A/a, which diverges for 
small necks. Thus, this spinodal cannot be reliably de- 
termined, if at all, by our numerical scheme. Within the 
approximation that the neck of the budded state shrinks 
to a size much smaller than the domain size, we find that 
buds are, in fact, metastable throughout the phase dia- 
gram. This is because, for small necks, there is a linear 
increase in the free energy due to the line tension when 
the neck expands, while both bending and surface tension 
contributions are quadratic. 

The remaining phase boundary shown in Fig. llOl is that 
of the transition from the flat phase (macro phase separa- 
tion) to the budded state. This phase boundary is simple 
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FIG. 10: The calculated phase diagram, indicating bud, 
caplet, and stripe phases, as well as the flat phase corre- 
sponding to macroscopic phase separation. This has been 
obtained by the numerical procedure outlined in Sec. [3 The 
solid lines indicate thermodynamic phase boundaries, while 
the dashed line indicates the limit of metastability of caplets. 
As discussed in the text, we expect that buds are metastable 
throughout the phase diagram, within the approximations 
used here. The thin region shown is that of the stripe phase. 
This vanishes at a triple point near (it, A) = (1,0.2). The 
bud, caplet, and flat phases meet at the other triple point 
along a — 0.5. Upper and right hand axes are dimensionless, 
while lower and left hand axes are in physical units for typi- 
cal values of the bending energy and spontaneous curvature, 
k = 100 pN nm,Co 1 = 100 nm. 

to estimate. The free energy per unit area of the bud is 
simply — ^kcq, while the bud accounts for an excess area 

of 4-7T (2/co) 2 , assuming an optimal mean curvature of cq. 
As noted above, the radius of the neck tends to zero. 
This results in a phase boundary given by a — 1/2 in 
reduced units. This simple phase boundary estimate is 
well borne out by the numerical minimization, as shown 
in the figure. 

VI. DISCUSSION 

The transition from caplets to stripes with increasing 
<f> and t can be understood in simple physical terms. On 
the one hand, with increasing <f> the unit cells in the stripe 
and caplet morphologies are of size ~ l/<f) and ~ \j\f<\>, 
respectively. Thus, with increasing (j), the constraints on 
the membrane profiles due to the unit cell boundaries 
become more costly in the caplet phase than the stripe 
phase. On the other hand, with increasing tension r, the 
regions of largest slope (i.e., largest difference between 
projected area and membrane area), which occur near 
the domain boundaries, become energetically more costly 
because of the work required against tension. For caplets, 
this represents a larger fraction of the domain area than 
for stripes. Thus, stripes are preferred for sufficiently 
large r. 

From the numerical approach of the previous section, 
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we see, perhaps surprisingly, that the Monge approxi- 
mation works rather well to describe the various phases 
and transitions. Specifically, we see from Fig. OH that for 
A < 0.2 there is little disagreement between the Monge 
approximation and the numerical solutions that can ac- 
count for finite slopes. This is particularly true of the 
caplet phase. 

In general, we expect the Monge approximation to be 
best for large a, which tends to flatten the membrane. 
Since the phase boundaries between caplets, stripes, and 
the flat phase tend to occur for a fixed r = a 1 / 2 A (specif- 
ically r ~ 0.21 and t = 0.25 for the caplet-stripe and 
stripe-flat boundaries, respectively, in the dilute limit), 
we expect the Monge approximation to be valid near 
both transitions, except for large A, which is what we 
find. Here, of course, a tends to be smaller. 

As noted above, however, the discrepancy is smaller for 
the caplet phase. The deviation that we find for stripes 
suggests an instability toward larger domains and, ulti- 
mately, the flat phase. This tends to be compensated in 
the case of caplets by the line tension A, which prevents 
domains from growing. This enhanced stability of caplets 
is also consistent with the observation from the phase di- 
agram in Fig. 1101 that the caplet-stripe phase boundary 
shifts in favor of caplets as A increases. We find, in fact, 
a triple point, where caplet, stripe and flat phases meet, 
although this is hard to see in the figure. For larger A, 
no stripe phase is observed. 

Away from these transitions, however, we begin to 
see significant deviation from the Monge results for the 
largest A as sigma decreases. Thus, the Monge approxi- 
mation tends to fail to the left in Fig. EH In particular, 
the Monge approximation is insufficient to characterize 
the transition between buds and caplets. 

This work shares some common features with sev- 
eral previous theoretical studies. Undulating stripe and 
hexagonal phases were predicted in the limit of weakly 
segregated two-component films based on a phe- 
nomenological coupling of local membrane composition 
and shape. These phases are analogous to the stripe and 
caplet phases discussed in this letter for well-segregated 
mixtures. The approach of Ref. [2^] has recently been 
utilized to study one-dimensional shape profiles of two- 
component membranes and vesicles in the strong segre- 
gation regime as well [2^,|2i|. The stripe phases of planar 
membranes in the limit of small domains or large bend- 
ing energies reported in Ref. 26] are qualitatively similar 
to our stripe mesophases. The stable stripe and caplet 
mesophases reported in this letter should occur in phase- 
separated two-component films under tension, provided 
that one phase has finite spontaneous curvature and that 
the line and surface tensions are sufficiently small. In 
particular, one must have r < 1/4 and A < a to observe 
these mesophases. The calculation presented above as- 
sumes an asymptotically flat matrix phase, as in the case 
of a globally flat monolayer or bilayer membrane. How- 
ever, our results should be qualitatively valid for two- 
component bilayer vesicles, provided that the vesicle di- 



ameter D is sufficiently large compared with the charac- 
teristic domain dimension, i.e., if CqD 1. 

Gozdz and Gompper performed a similar calculation to 
ours in a different "slice" of parameter space. They con- 
sidered phases with different signs of spontaneous curva- 
ture, and predicted a variety of different budding, stripe, 
and circular domains. It is difficult to compare the re- 
sults directly, since they chose only one parameter set, 
A = 0.5, a = 0.4. For our model, with just one spon- 
taneous curvature, this leads to a budded state. They 
show that this set of parameters can lead to a variety of 
phases that can be stabilized by the different spontaneous 
curvatures. 

We estimate the region of domain stability for the fol- 
lowing values of the parameters in our model: 

K = 25k B T ~ lOOpN nm, c 1 = lOOnm. (6.1) 

The results are shown schematically in Fig. QJJ Thus, 
caplet domains of order 100 nm in size (as in Ref. |40|) are 
stable for line tensions A < 0.3 pN and surface tensions 
a ~ 10~ 2 mN/m. For comparison, membrane rupture 
occurs for surface tensions of order lOmN/m 

Finally, we wish to note a couple of important gen- 
eral observation relevant to phase separating membrane 
systems. First, given the coupling between membrane 
composition and shape, it is possible to have both 
composition- induced changes in shape or curvature pH 
l29l l31l. as well as curvature-induced changes in composi- 
tion [3f| . However, unless the curvature is of a molecular 
scale (i.e., if curvature radii are comparable to the size a 
of a single lipid), it is generally expected that composi- 
tion will drive shape, or, in other words, that chemistry 
determines geometry. This is because, on a per molecule 
basis, the curvature energy ko 2 Cq is very small compared 
with kT . Thus, unless one happens to be very close to 
a critical point |30|. the shift in chemical potential due 
to the local curvature of the membrane is expected to 
be insufficient to cause a significant composition change. 
Second, the Gaussian curvature and associated modulus 
K, which we have so far neglected, is expected to affect 
the stability of the caplet and budded phases described 
above. In fact, this can have an order one effect on both 
phases, since the corresponding shift on the curvature en- 
ergy is of order Akc 2 , per unit area, where AR represents 
the difference in moduli between the two phases. This is 
comparable to the (mean) curvature energy contribution 
in both phases. There will be no effect, however on the 
stripe phase free energy. Thus, the main qualitative effect 
of taking into account the Gaussian curvature effects will 
be on the stability of the bud and caplet phases relative to 
the stripe phase. For instance, the expected increase or 
decrease in the energy of bud and caplet phases may re- 
sult in a more pronounced or even absent stripe. Finally, 
we have imposed the strong segregation limit, in which 
the compositions of the two phases are unaffected by the 
local bending. This is likely to break down near critical 
points, at which point the compositions will vary in con- 
cert with the degree of bending 

ins Hummus. 
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In this case the local line tension can be expected to be 
lower and depend on the degree of local membrane defor- 
mation. Moreover, elastic constants and the spontaneous 
curvature will be inhomogeneous in the more highly de- 
formed regions. Possible effects of this could be to stabi- 
lize the caplet phase in favor of the stripe phase (since line 
tension may be reduced), as well as shifting other phase 
boundaries. However, as we have emphasized above, the 
strong segregation limit is expected to apply more gen- 
erally, because of the wide separation of scales between 
the molecular and curvature lengths. 
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APPENDIX A: DERIVATION OF FREE ENERGY 
IN MONGE GAUGE 

In this appendix we derive the free energy and bound- 
ary conditions for the Monge gauge calculation. 



1. Free Energy Variation 

Consider an impurity phase arranged into domains. 
The free energy is given by Eq. (|2.1|l . The bending free 
energy with respect to a flat reference state is given by 



Ft, 



dS 



c(r)-2c(r)c + 4 l-ft dS c(rY , 



Sa 



where c(r) is the mean curvature, Co is the spontaneous 
curvature of the impurity phase, k is the bending modu- 
lus, and dS is the surface area. The structure has a frame 
tension a against which the membrane does work, with 
energy 



F s = a{A fr - f dS x \, 

I JSa+Sb J 



(A2) 



where dS± is the projected membrane area and Af r is 
the frame area. The line tension energy is proportional 
to the interface between phases A and B and is given by 
Eq. |£3J>. 

The membrane composition <p is given by 



Sb + Sa ^ 
The two areas Sa and Sb are determined by the frame 

tension and a Lagrange multiplier fi to control the im- 
purity phase B. Since we assume a constant area per 
lipid, (i is proportional to the chemical potential for B, 
or more strictly the exchange chemical potential differ- 
ence between A and B. We must minimize 



G (M,M) = F b + F s + Ft + ti [ dS 

Js B 



(A4) 



over variations in the surface shape, h(r) — > h(r) + 5h(r), 
where the line tension energy Fi is given by Eq. (|2.3|l . 
We will later minimize the resulting free energy over the 
domain dimensions. 

In the Monge gauge in the 2d geometry of the caplet 
case, 



1 d ( dh 



dS 
dS ± 



l + \(Vhf 



1/2 



d 2 r 



d\ 



(A5a) 

(A5b) 
(A5c) 



(Al) and the variation of G over h leads, to lowest order in h, 



5G 
2^ 



Sh(r ) f hr (r ) + Sh(R) f hR (R) + Sh'(r ) r r (r„) + 8h'(R)T R (R) + / d 2 r6hff^ + 



d 2 rShf^_ 



(A6) 



in terms of the following forces / and torques T: 

f hr = 2n r [k (c' a - c B ) - (a + p,)h' B + ah' A ] (A7a) 

f hR = 2tt R[ah' A - kc' a ] (A7b) 

r r = 2tt nr (cb — ca — Co) (A7c) 

T b , = 2ttRkc a . (A7d) 



Here we have defined ft, = fi + kc 2 , which is consistent 
with a complete expansion to second order in h |3?t| . The 
Euler-Lagrange equations, which determine the profile 
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h(r) are 

Je-l 
Je-l 



k(c' b 
k(c' a 



rc'i) 



(a + fl)rc B 
arc a = 0. 



(A8a) 
(A8b) 



Because the free energy and the forces and torques de- 
pend only on derivatives of h, we will solve for the slope 
i] = dh/dr instead of the height profile. 

Eauations lA7IA8l are specific for the cylindrical geome- 
try of the caplet morphology. The analogous set of equa- 
tions for the stripe morphology is: 

fhr 

L 

fhR 

L 

El 

L 
Tr 
L 

,B 

Je-l 
,A 

Je-l 

where now h is a function of the cartesian coordinate x, 
and the curvature is given by c = d 2 h/dx 2 . 



2. Boundary Conditions 

The boundary conditions are specified by the forces 
and torques applied to the membrane. There is no ap- 
plied torque at the interface, rp, and we assume no ap- 
plied torque at the boundary j42|, so 
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(A9e) 


KC'I- 


ac A 


= o, 




(A9f) 



r r (r ) = 
T R (R) = 0. 



(AlOa) 
(AlOb) 



In this Appendix we will explicitly describe the boundary 
conditions for the two dimensional caplet geometry; the 
same conditions hold for the one dimensional stripe ge- 
ometry, with r — ► xq, unless noted. The second condition 
specifies zero curvature on the boundary while the first, 
first given by Kozlov and Helfrich, relates the difference 
in the mean curvature to the spontaneous curvature of 
the inner surface. Next, we assume there are no vertical 
forces applied to the membrane, 



fhr(ro) 

fhR(Ro) 



(Alia) 
(Allb) 



However, it is straightforward to show (for both the one 
and two dimensional cases) that these vertical bound- 
ary conditions are in fact proportional to first integrals 
of the Euler-Lagrange equations. Specifically, in two di- 
mensions, 



fh(r) = r (kc' (r) - ah! (r)) 
"drf e -i[h(r)) 



constant. 



(Al2a) 
(A12b) 
(A12c) 



Since fh(r) is constant for all r, the boundary condition 
at the junction, Eq. l|Alla() . is irrelevant and we are left 
with Eq. l|Allb|) . A similar result applies for the one 
dimensional stripe geometry. 

To completely specify the problem we need a few more 
conditions. Symmetry about the domain center requires 



h B (0) = h"'{0) = 0. 



(A13) 



Another condition on the derivative of the profile h'(ro) 
results by requiring a well-defined energy. If hi has a 
discontinuity at the boundary, then the curvature has a 
delta function as the boundary is crossed. This leads to 
a non-physical and singular energy in the region of the 
junction r = r$. Hence we are led to our final boundary 
condition, 



ti B ( r o) = h' A {r )- 



(A14) 



For completeness, we reproduce the relevant bound- 
ary conditions that remain after recognizing the redun- 
dance of Eq. (|Allap . For the single-domain problem, 
from Eqs. (|A7blA7d|) we have: 



CB - CA 



r/i A (oo) — kc a (oo) 
ca(oo) 



co 



0. 



(A15a) 
(A15b) 
(A15c) 



The three equations above, together with the three 
conditions of Eqs. I|A13IA14|I . give us six conditions. The 
Euler-Lagrange equations (Eq. IA8|) are second order dif- 
ferential equations for c B (r) and ca(^), which have four 
constants. The solutions to these equations must be in- 
tegrated once to obtain rj(r), which yields six constants. 
Finally, we determine the domain size ro by minimizing 
the total free energy (after substituting the solution into 
the original free energy) over ro. Hence we have a closed 
set of equations which determine the height profiles h B (r) 
and hA^r) and the impurity domain size ro, for a given 
(arbitrary) frame area. 

For the multi-domain problem the boundary condi- 
tions must be modified. We still have Eqs l|A13IA14() 
above, but instead of boundary conditions at R = oo, 
we impose a symmetric profile at the Wigner-Seitz cell 
boundaries: 



4(i?) = h' A (R) 
replacing Eqs. HA15blA15c(l . 



0, 



(A16) 



APPENDIX B: SOLUTIONS 
1. Stripes 

For stripe domains the Wigner-Seitz cell is an impurity 
stripe of width 2xq embedded inside a majority stripe 
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with total with 2Xo. The curvature is c = d 2 r]/dx 2 , and 
the Euler-Lagrange equations, Eqs. (|A9() . are 



d 4 h A d 2 h A 

K — : — -. a , » = 0, 



dx 4 

{(T + jl)'- 



dx 2 



(Bla) 



0. 



dx A dx 2 
To calculate the free energy we only need the slope, 

dh 



(Bib) 



V 



dx 



(B2) 



Eqs. (|B1J) may be integrated once, with integration con- 
stant zero (all odd derivatives d^h/dx n vanish at the 
origin). This leaves the following Euler-Lagrange equa- 
tions: 

d 2 "q A 



d(, 2 



077A = 0, 



d 2 m 
dc 2 



(<t + p)rjB = 0, 



(B3a) 



(B3b) 



where we have introduced the following dimensionless 
variables: 



(B4a) 
(B4b) 

(B4c) 
(B4d) 

Note that /2 and a both have dimensions of energy/area, 
because we have implicitly assumed an incompressible 
membrane with an unchanging area per molecule. The 
solutions are 

sinhCsC 



XC 
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2 



2 

KC 
XqCq. 
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sinhCA(S - Co) ' 



where 



no 

ci 



[Cb coth CbCo + Ca coth Ca(S - Co) 

U + jl 



(B5a) 
(B5b) 

(B6) 
(B7) 
(B8) 

The hyperbolic functions satisfy the conditions of zero 
slope at the domain center C = and domain boundary 
£ = S (Eq. HA13I IA16j) . while the amplitude % is de- 
termined by the discontinuity of curvatures at £ = Coj 
according to Eq. (|A15a(l . 

After some rearrangement, the free energy (Eq. IA4|) 
may be written as 

2L^r A+ ^°-y A (B9) 

The number of stripes is given by TV = ^Jj; an d the total 
free energy is G — NGq. We can thus write the free 
energy per frame area <?, in dimensionless form, as 
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G 



1 



Af r £1 



A + pCo 



»7o 
2 



(BIO) 



2. Caplets 

We consider circular domains of inner radius tq and 
outer radius R. The Euler-Lagrange equations for the 
circular geometry, Eqs. (|A8(I are given by 



n{c A + rc A ) — arc A = 0, 



(Bll) 



and similarly for cb- In cylindrical coordinates, the mean 
curvature is 



_ ld_ (dh 
r dr V <9r 



(B12) 



Equation QBllfl may be integrated once to give 

r {kc' — ah'} = constant. (B13) 

Again, h'(0) — c'(0) = by symmetry, so the constant 
vanishes. As with the stripe case, we only need the slope 



dh 

in terms of which Eqs. (|A8(I become 



(B14) 



d 2 r lz 1 % / 2 1 \ 

for i = A, B, where p = tcq. Implementing the boundary 
conditions yields the following profiles: 



Vo Ii[(bPo] 

Vo K i [Capo] -Ah [CaPo] ' 
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Ca/i[CbPoJ k i[C aH-^/i[Capo] + CbJoICbA) 



(B16a) 
(B16b) 

(B17) 
-B18) 



where i? = i?co and h, h, Ko, and -ftTi are modified Bessel 
functions . The free energy per caplet is given by 

G 1 

ir— = -Vo + {PP - Vo) ■ (B19) 

The number of caplets is N = Af r /(irR 2 ), so the total 
free energy is 

° = ^ r = w [2X+p ^- Vo] - (B20) 

Eqs. (|B20I IB10|) are used in Eqs. (|2.8|l to evaluate the 
equilibrium domain sizes. 
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